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TECHNICAL MEMORANDUM X-64921 


A TIMELINE ALGORITHM FOR ASTRONOMY MISSIONS 

I. SUMMARY 


An algorithm is presented for generating viewing timelines for orbital 
astronomy missions of the pointing (nonsurvey/scan) type. The algorithm 
establishes a target sequence from a list of candidate targets in a way which 
maximizes total viewing time. Due to lighting constraints, dim targets are 
scheduled only during the antisolar portion of each orbit. They require long 
observation times extending over several orbits. The viewing of one dim target 
is completed before the next dim target is scheduled. Brighter targets, which 
have less restrictive lighting conditions, are scheduled during the portion of 
each orbit when the dim target cannot be viewed. Usually, the viewing require- 
ments of two or more bright targets can be completed during this part of the 
orbit. An enumeration algorithm described in the appendix is employed to 
select these targets. A minimum slew heuristic is employed to select the 
sequence of dim targets to be viewed. A complete description of the algorithm 
is provided by Section II. 

Application of the algorithm to the Large Space Telescope (LST) pro- 
gram is discussed in Section III. The LST timelines illustrate the full capability 
of the algorithm and how it can be used to simulate and/or schedule complex 
space astronomy missions. 


II. INTRODUCTION 


The timeline algorithm presented in this report was originally developed 
for the LST program but is now applied also to other spacecraft and astronomy 
missions of the pointing (nonsurvey/scan) type. This algorithm determines a 
target sequence, or viewing path, from a list of candidate sources in such a 
way that total viewing time over a specified portion of the mission is maximized 
while maintaining a uniform distribution of target magnitude. 



The selection of an optimal path is typically restricted by orbital lighting 
conditions, by the time intervals each star can be seen, by the observation time 
required for a given star and sensor, and by observational philosophy itself. 

For example, viewing requirements for faint targets are normally longer than 
one interval, so a number of revolutions are required to complete observation. 
In addition, stray light constraints make it preferable, if not necessary, to view 
dim sources on the antisolar side of an orbit. This may be either strictly the 
orbital shadow interval or, in the case of telescopes equipped with sunshades 
and baffling systems, a longer interval centered around orbital midnight. 
Philosophy enters at this point in formulating approaches to the space observa- 
tory scheduling problem which will yield high mission utility, low program cost, 
and maximum scientific return. Subsystem capabilities, instrument sensitivi- 
ties, and target dwell times are the primary factors that influence the selection 
of a specific approach. 

The LST, for instance, is designed to observe very dim sources, which, 
in turn, requires very long observation times extending over many revolutions. 
The large dynamic range which makes this possible, however, results in very 
short times for relatively bright objects. The problem of how to mix such a 
wide diversity of times in the most eiiicient manner is the key to the logic 
developed in the algorithm. As originally conceived, the LST had the ability 
for relatively fast slewing, making it possible to supplement the viewing pro- 
gram with brighter targets during periods of prime region occultation. This 
allowed a separation of the brighter targets from the very faint into two distinct 
portions of each revolution; i. e. , the portion centered around orbital midnight 
is reserved for faint sources and afl other sources are scheduled within the 
remaining portion. Restricting observations in this manner is equivalent to 
dividing the very long viewing times from the short. 

Since the viewing of faint sources is one of the primary LST objectives, 
any scheduling technique should provide as much faint source viewing as 
possible. This can be accomplished by beginning observation of a dim target 
as soon as it is viewable each dark-side pass and continuing observation of the 
current dim target until it is lost from view. In addition, by applying the time 
restriction from the previous paragraph and by slewing to subsequent dim 
targets requiring the minimum maneuver time, primary viewing will attain the 
highest possible efficiency. Since slewing between primary targets would then 
be infrequent, this portion of each revolution can be made very efficient. 

This scheme implies, of course, that either there is little interest in 
observing brighter sources or that the spacecraft/instrument platform has a 
sufficiently high slew rate capability to provide secondary viewing on the bright 
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side of each revolution. A number of Spacelab payloads and the 3 m class LST 
design do have this capability. The current LST redesigns with their small and 
inexpensive pointing and control systems, on the other hand, render secondary 
viewing all but impossible. In this case a somewhat different approach must 
be taken. 

Although the scheduling algorithm can be applied more generally, the 
specific approach modeled in this report follows the secondary viewing scheme 
described above. Target acquisition/loss intervals based on stray light con- 
straints and occultations are calculated by another routine and fed into the 
scheduling model by tape. Candidate targets are divided into two classes for 
each instrument complement: one for faint sources viewed on the antisolar side 
and one for supplementary bright-side viewing. The telescope returns to the 
current primary target every dark -side pass until the required observation 
time is satisfied and then slews to the next dim target requiring the smallest 
maneuver time. Observation times are calculated bj^ summing integration times 
from input curves over the instrument complement scheduled. Secondary view- 
ing, on the bright side, is determined by the algorithm, which maximizes total 
viewing time while minimizing losses due to slewing and unused time. The slew 
from the primary to the secondary region is not permitted until the primary 
target is lost from view, and the telescope must slew back and be ready to view 
the primary target again (or next dim target) as soon as it is available. The 
current primary target then sets the begin and end point of the bright-side 
target path. 

The algorithm, therefore, performs two distinct tasks: (l) selects a 
sequence of faint sources over a given mission segment and (2) determines 
optimum bright-side viewing paths that fit between the dark-side viewing paths. 
The latter is achieved by implicitly enumerating all possible and promising 
paths, making use of fathoming tests for eliminating nonpromising paths to find 
the one with the best objective value. Notice that the primary target path is 
independent of the bright-side path and is, therefore, the same regardless of 
the presence or absence of secondary viewing. In addition, bright side paths 
have no effect on each other (except the diminishing list of candidate targets as 
they are scheduled) , so substitution or deletion of scheduled targets does not 
require a completely new timeline. 

Although this scheme produces high observation efficiency, it does have 
several drawbacks. For example, the slews within a target region can be made 
quite small, but the two slews per revolution connecting the primary and 
secondary regions are of necessity relatively large to avoid earth occultation. 
This means that the pointing/control system must provide much more slewing 
(magnitude and frequency) than would be necessary for dark-side-only viewing. 
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In the case of free flying satellites, such as LST, this is a significant design 
factor but, for pallet -mounted telescopes on the Shuttle, this is of little conse- 
quence. Another disadvantage is some additional complication in data manage- 
ment, since the approach implies a separation of ,T dim" and "bright” target data. 
The impact of this problem, however, is softened if bright-side targets are 
restricted to single^pass observation, so only the faint source memory needs to 
keep track of subsequent revolutions. 

In cases where this "multitarget" viewing is not attainable (e. g, t the 2 m 
class LST redesign offers a slew rate so slow that supplementary viewing is 
virtually eliminated), the dark-side viewing interval should be "stretched" as 
much as possible with protection devices such as sunshades and baffles and 
efficient use should be made of the available time. The remainder of each 
revolution would then be used for data transmission, subsystem monitoring, 
solar array reorientation, antennae positioning, instrument preparation, etc., 
and any other spacecraft functions that are inhibited during observation. Since 
"short" and "long" viewing times are no longer divided into separate intervals, 
efficiency can be maintained by determining an acceptable range of "mixture 
ratios" (short observations /long observations) which will keep slew losses to a 

]j!51uIluUni * 


Other schedule philosophies are possible but will not be discussed here 
since they are beyond the scope of this report. 

The sequence of faint sources is determined by a simple heuristic. The 
first target in the sequence is either specified by the user, is the earliest 
available source, or is determined from the position of the previous target. 

The next target is selected by listing all the candidate sources that are visible 
at the time viewing of the first source will end and, from this set, choosing the 
target requiring the minimum slew time. In general, then, the (k + l) target 
in the sequence is the target requiring the minimum slew time from the kth 
source which also belongs to the set of targets in view at the time viewing of the 
kth source is completed. 

Since the faint target path is a relatively slow moving affair, the total 
efficiency of each revolution will depend primarily on the secondary viewing 
path selection. An optimization procedure has been developed to assure that 
the selected target sequences achieve the highest possible efficiency. The 
techniques of integer linear programming are used to search the feasible viewing 
paths to find the one which maximizes path value. An input value is assigned to 
each candidate source and is used in a function which determines path value. 

This function increases with weighted observation time and decreases with slew 
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time and waiting time (unused time at end of bright side) . Any weighted values 
can be assigned to each target, but the experience gained from a statistical 
analysis of this parameter with random target samples and a number of mission 
studies indicates that equating the weight to required viewing time results in 
both the maximum total viewing time and a uniform, unbiased distribution of 
source magnitudes (equivalent to observation times). The statistical analysis 
also showed that the function described above consistently out-performed and 
was less biased than functions which either maximized total viewing time or 
minimized slewing time only. 

A detailed description of the optimization technique is provided in the 
appendix. 


III. TYPICAL APPLICATION OF ALGORITHM 


One of the first applications of this algorithm was in timeline simulations 
for an LSI' Design Reference Mission (DRM) document which was to have been 
published in September 1974. Although the LST program experienced a major 
redefinition shortly before this time and the document was never published, the 
timelines generated for it illustrate the full capability of the algorithm and the 
manner in which it can be used to simulate and/or schedule complex space 
astronomy missions. The input conditions, instrument and target modeling, 
schedule rationale, sample timelines, and some of the statistical results are 
discussed in this section. 


A. Input Assumptions and Instrument Characteristics 

The first step in the simulation process was to set down the operational 
assumptions and construct the instrument package math model used by the 
computer. A total of nine instruments were defined, eight of which were 
science instruments, while one, the target acquisition camera (Slit Jaw Camera), 
was an auxiliary element. Characteristics of these fictitious instruments are 
presented in Table 1. The arbitrary devision of the spectrographs into pairs 
was motivated by a desire to simulate sensors with different spectral ranges. 
Estimated integration time curves for each instrument are shown in Figures 1 
through 4. These curves, which are input to the program by curve fit equations, 
define how long given sensors must observe targets of given magnitude to obtain 
one high quality data frame. Times ranged from a maximum bf 10 hr to < 1 min. 
A constant 10 min operation was assigned to the Slit Jaw Canfera (SJC), which 
was required with any complement including a spectrograph. 
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TABLE 1. INSTRUMENT CHARACTERISTICS ASSUMED 


Instrument Name 

Abbreviation 

Magnitude 

Range 

Power 

(watts) 

Maximum 
Data Requirement 
(bits /fra me) 

High Resolution 
Camera (f/24) 

HRC 

^ 25 
25-30 

90 

105 

2 x io 8 

Faint Object 

FOS-1 

10-19 

105 

2 x lo 8 

Spectrograph 

FOS-2 

12-24 

135 

2 x l() 8 

High Resolution 

IIRS-1 

7-14 

105 

2 x I f/ 

Spectrograph 

HRS -2 

7-16 

135 

2 x io 8 

Photometer 

PHOT. 

S 30 

28 

6 x 10 7 

Astrometer 

ASTR. 

^ 20 

15 

5 X 10 7 

| Infrared 

IK 

% 22 

2.2 

1 x io 1 

Slit Jaw 

Camera (Target 
1 A equisition 
| Camera) 

SJC 

(TAG) 

£ 30 

40 

2 x HP 


It was assumed that the high resolution camera (HRC) could operate 
simultaneously with any other science instrument but that all other sensors in 
a specific complement would operate in series. In other words, associated with 
each target is at least one ai’ray of sensors, each of which will record a data 
frame. If the HRC is involved, it will operate at the same time with one of the 
other instruments. Thereafter, the remaining sensors will operate one at a 
time until all desired data are collected on the target. In scheduling observa- 
tions, the program will sum the integration times for each instrument in the 
complement, adding any specified time gaps between instruments (for switch- 
ing, calibration, filter change, etc. ) . The total time resulting from this calcu- 
lation is the time that the scheduler routine must allow for this target. With 8 
science instruments there are 255 possible complements, the number of 
instruments in a complement ranging from 1 to 8. Of this total, however, many 
seem unlikely and would be rarely, if ever, used. For this reason a total of 33 
complements with up to 5 instruments each were selected; these complements 
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Figure 4. Integration time curves — IR instrument 
photometer, and astrometer, 
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are defined in Table 2. This list is by no means exhaustive but covers a wide 
range of possibilities. Each grouping appeared likely for relatively frequent 
use. 

TABLE 2. INSTRUMENT COMPLEMENT MATRIX 
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B. Target Modeling 

The 350 candidate targets fed to the computer were all pseudo-targets 
generated by an auxiliary random number program; i.e. , any resemblance 
between one of these targets and a real heavenly body is purely coincidental. 

Use of pseudo-targets has a number of advantages: (l) any model of the celes- 
tial sphere or portion thereof can be simulated without bias, (2) it is much 
faster than digging from star catalogs, (3) large gaps in existing data can be 
filled (e.g. , very faint objects), and (4) it eliminates the emotional involve- 
ment sometimes encountered with real heavenly bodies. 

Fifty of the targets were designated as "faint” and generated in the anti- 
solar region depicted in Figure 5. The remaining targets were ’’bright” and 
taken from the rectangular areas close to the sun constraint circle also shown 
in the same figure. These three sky areas represent the best ’’viewing windows” 
for faint and bright targets, respectively, during the 2 day period of the 
enclosed timelines. 

All objects were invented at random locations within each area and given 
random integer visual magnitudes. ’’Brights” were assumed to vary between 
m = 10 and m = 22 , while ’’faints” varied from m =23 to m =29. No 

v V y V 

magnitude brighter than 10 was used because (l) integration times at this mag- 
nitude for most of the sensors are already g 1 min and (2) the LST will prob- 
ably have little interest in objects brighter than this. One-third of the pseudo- 
targets were randomly classified as UV only* one-third as IR only, and one- 
third observable throughout the LST spectral range. 


The final step in the input construction was to assign an instrument 
complement to each target. This was done arbitrarily* but within spectral con- 
straints so that each complement was represented by at least one target. Each 
of the targets scheduled in the timelines is defined in Table 3. 


C. Operational Mode Scenario and Rationale 

The isochronal viewing contours shown in Figure 6 define the total time 
per orbital revolution that an object can be viewed within stray light constraints. 
These constraints are determined by the amount of extraneous light that can be 
tolerated at the focal plane and maintain and image of given quality for a specific 
sensor. The amount of stray light at the focal plane, in turn, is a function of 
star magnitude, baffle design, sun shade design, zodiacal light level, and 
vehicle orientation relative to the sun, moon, and bright earth. 
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Figure 5. General regions of interest on the celestial sphere and faint and bright target areas used in timeline. 









TABLE 3. PSEUDO-TARGETS SCHEDULED, OPTIONS 1 AND 2 
(Listed in Scheduled Order) 


Target No. 

R.A. (deg) 

Declination 

(deg) 

Visual 

Magnitude 

Spectral 

Type 

Instrument 

Complement 

50 

62.092 

77.526 

28 

IR 

Q 

130 

221.533 

29.021 

11 

uv 

A 

263 

317.608 

4. 151 

13 

Comb. 

Y 

113 

224.235 

27.948 

13 

IR 

Y 

209 

318.520 

-2.511 

15 

UV 

Q 

237 

317.357 

0.257 

11 

IR 

E 

66 

213.839 

22.756 

15 

UV 

Q 

72 

220.797 

23.147 

11 

UV 

U 

305 

331.809 

2.678 

T5 

Comb* 

X 

59 

220.539 

22. 116 

15 

uv 

Q 

71 

221.218 

23.823 

11 

Comb. 

X 

216 

324. 570 

-4. 016 

15 

UV 

Q 

155 

208.280 

31.282 

15 

IR 

X 

214 

322.772 

-4.134 

15 

UV 

F 

215 

323.577 

-3.833 

22 

UV 

G 

64 

1 

215.634 

22.065 

15 

uv 

F 

226 

325.089 

-3. 296 

14 

uv 

Q 

272 

328.203 

1.331 

15 

IR 

E 

73 

219.668 

23.927 

15 

IR 

T 

319 

323.905 

6.162 

12 

Comb. 

Y 

3 

77.760 

56.802 

24 

IR 

Q 

169 

218.361 

32.136 

12 

Comb. 

Y 

321 

326.428 

5.781 

14 

UV 

Q 

76 

213,551 

23.602 

14 

UV 

Q 

247 

322,841 

-0. 966 

12 

UV 

Q 


14 










TABLE 3. (Continued) 


Target No. 

R.A. (deg) 

Declination 

(deg) 

Visual 

Magnitude 

■ 

Instrument 

Complement 

298 

324.013 

5.697 

18 

uv 

G 

12 

78. 189 

58.289 

24 

IR 

Q 

56 

195.464 

21.039 

12 

Comb. 

Y 

304 

330.007 

3.001 

11 

Comb. 

Y 

111 

197.603 

26.638 

14 

UV 

X 

244 

316.972 

1.425 

12 

UV 

U 

224 

321.331 

-1, 812 

17 

Comb. 

S 

13 

76.145 

58.402 

26 

IR 

Q 

171 

215.413 

32.629 

14 

Comb. 

Q 

225 

322.583 

-3.071 

17 

Comb. 

s 

212 

321.644 

-3.329 

16 

IR 

E 

154 

210.406 

31.089 

11 

Comb. 

Y 

230 

335.863 

-6.204 

13 

UV 

F 

97 

216.714 

25.076 

12 

UV 

Q 

62 

' 217.071 

22.448 

16 

IR 

E 

337 

317.357 

1 

0.257 

11 

IR 

T 

99 

213.990 

25.812 

16 

Comb. 

Q 

222 

318.673 

-1.464 

11 

UV 

U 

221 

318.562 

-0.694 

11 

Comb. 

X 

21 

73.313 

60.757 

23 

IR 

J 

187 

196.220 

33.792 

12 

Comb, 

T 

261 

339.848 

-3.967 

14 

Comb. 

X 

80 

201.407 

23.654 

13 

Comb. 

F 

283 

335.919 

-0. 504 

12 

Comb. 

X 

284 

337.636 

-0. 813 

15 

IR 

E 


15 








TABLE 3. (Continued) 


Target No. 

R.A. (deg) 

Declination 

(deg) 

Visual 

Magnitude 

Spectral 

Type 

Instrument 

Complement 

94 

223.713 

25.188 

12 

IR 

U 

206 

340.807 

-9.821 

12 

Comb. 

Y 

74 

217.955 

23.766 

17 

Comb. 

S 

89 

218.696 

24.260 

15 

IR 

E 

286 

341.935 

-1.888 

12 

UV 

Q 

75 

216.208 

23.006 

17 

Comb. 

s 

95 

221.736 

25.281 

15 

IR 

E 

307 

I 

332.962 

2.273 

11 

UV 

U 

275 

334.012 

-0.485 

20 

UV 

G 

163 

1 on 

J- M V t 

Ol AOO 
tJ X • WIJ ^ 

i /\ 
111 

Cornu, 


342 

331.479 

5.757 

11 

UV 

U 

343 

333.088 

5.384 

14 

IR 

E 

133 

202.754 

29. 225 

12 

Comb. 

X 

234 

344.080 

-7 . 588 

11 

UV 

U 

233 

342.403 

-7 . 143 

13 

IR 

E 

20 

74.164 

60.375 

26 

IR 

Q 

136 

195.690 

29.048 

12 

UV 

Q 

313 

346.770 

-0.519 

10 

UV 

Q 

192 

L 

209.674 

34.148 

11 

UV 

U 

165 

222. 116 

32. 951 

10 

IR 

R 

231 

336.237 

-6.731 

11 

Comb. 

X 

i 

j 232 

339.407 

-7.103 

19 

UV 

G 

i 

157 

206.855 

31.199 

11 

UV 

U 

227 

329. 583 

-4.744 

17 

Comb. 

S 

217 

327.997 

-5.545 

14 

IR 

E 
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TABLE 3. (Concluded) 


Target No. 

R.A. (deg) 

Declination 
( deg) 

Visual 

Magnitude 

Spectral 

Type 

Instrument 

Complement 

84 

192.370 

23.735 

11 

UV 

U 

220 

345.217 

-8.990 

17 

Comb. 

S 

19 

77.054 

60.845 

26 

IR 

Q 

81 

200.888 

23.219 

11 

UV 

X 

79 

203.010 

23.881 

14 

IR 

E 

315 

321.233 

8.073 

10 

UV 

R 

254 
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The exact size, shape, and position of the isochronal contours varies, 
but they all show that the maximum available viewing time per revolution for any 
object exists in a region centered around the antisolar point and that the maxi- 
mum time an object can be continuously observed per revolution is approximately 
45 to 52 min (depending on magnitude). The antisolar point, of cours.e, moves 
across the sky at the same rate as the sun and is opposite in declination and 
right ascension. The size of the maximum viewing region decreases as object 
magnitude becomes fainter, and contours outside this region tend to merge with 
it. 

Also moving with the sun is the zodiacal light, the primary source of 
unwanted diffuse background light [ l] . At any given time two symmetric 
regions exist centered <^65 deg north and south of the antisolar point, where the 
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Figure 6. Aitoff -Hammer projection of isochronic viewing contours in ecliptic 
coordinates for 25th magnitude targets on day 1. 



zodiacal light levels are at a minimum. When possible, viewing of faint objects 

outside the two S , . . = 100 contours should be avoided, since this is the 
10(vis) 

approximate level at which integration times begin to increase as a result of 
background light. These regions are roughly circular, subtending half cone 
angles of about 40 deg and are symmetric with respect to the ecliptic plane. 


If these contours are superimposed on the isochronal line, the maximum 
viewing region is cut into two smaller regions north and south of the antisolar 
point. Prime time viewing of faint objects, then, is most efficiently accom- 
plished on the antisolar side of the orbit and within one of these resultant regions. 
This is especially important since an increase in the already long integration 
times required for faint sources will accumulate a large amount of •wasted time 
and a possible degradation in the data. For example, a 29th magnitude star 
viewed at v = 150 may require an estimated 15 to 20 percent longer 

observation time. 


The natural progression of the sun will favor different sections of the sky 
at different times of the year. Targets outside of the area swept out during the 
yearly cycle, if any, would be viewed when they come closest to one of the prime 
viewing areas. 


D. Sample Timelines 

Sample pages from two LST timelines are presented in Figures 7 and 8 
to illustrate program capability and the type of information that can be derived 
from a timeline analysis. Both timelines represent the same mission segment 
and the same input conditions. The only difference between them lies in the 
number of secondary viewing regions employed. 

Option 1, dubbed TL-4 (timeline #4) , scheduled targets in the prime 
(A) and secondary (b) areas depicted in Figure 5. Option 2, or TL-5, made 
use of a third area (c) on the opposite side of the sun from area B. Both 
secondary regions are "twin brothers" and are symmetric with respect to the 
orbit plane. To eliminate any bias in scheduling between these two areas, the 
pseudo-targets generated for B were mapped into C and assigned the same 
instrument complements. 

The reason for option 2 was that the stray light constraint characteristics 
that were modeled did not permit target intervals throughout the entire bright 
side of each revolution. Region B will provide observation opportunities through 
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only the first three-fifths of the bright side; similarly, region C covers the last 
three-fifths, so that the use of B and C together gives continuous coverage for 
secondary viewing. Since these areas required a major connecting slew, the 
objective was to compare the differences in efficiency and slew time between 
the two timelines. 


Each timeline is presented in a three page format. The first page con- 
tains the target schedule, orbital reference data, radiation passage, and 
instrument operations; the second page depicts spacecraft attitude, slewing 
information, and a revolution-by-revolution observation efficiency histogram; 
the final page shows the power and data rate profiles resulting from the 
instrument operations. 


Orbital revolution number shown on the first page is referenced from one 
ascending node to the next. The numbers on the target schedule line specify the 
target under observation. For the characteristics of these pseudo-targets refer 
to Table 3. The instrument complement codes in the tables are defined in the 
matrix of Table 2. The fourth line depicts the prime viewing interval of each 
■revolution; the solid lines enclose the antisolar side of th 


U .1 JLl- „ -1 - 1 
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lines bound the actual shadow interval. Line five indicates when the LST passes 
through the 10 protons /cm 2 /sec, = 50 MeV level of the South Atlantic Anomaly. 
The last group of linos on the page show which instrument(s) is taking data. 


LST interial attitude is depicted on the second page in terms of right 
ascension and declination of the target being observed. The next two lines give 
the size and duration of each slew maneuver and a running average of slew 
angle. All slews were cigenaxis since this type of maneuver requires minimum 
time and momentum. A 2 min settle time was deducted after every slew to 
allow for fine alignment and damping of transient motion. The final line on this 
page is a revolution-by-revolution histogram of observation efficiency; i.e. , 
the percentage of the orbital period spent observing. Each revolution here is 
referenced from the beginning of one antisolar interval to the next. 


The last page of the timeline format depicts the power and data rate pro- 
files resulting from the observatory schedules. A standby power of 1500 watts 
was assumed for the total LST system, so that the variation in the profile is 
due to the observation program. Since all instruments require long warmup 
times and since their standby power is low, all instruments were assumed to be 
on during the entire mission segment. The data profiles were based on the total 
bits required for a frame of data on each instrument given in the last column of 
Table 1. These numbers are divided by the observation times to determine an 
average rate that will produce the desired number of bits at the end of the 
observation. Other buildup schemes could be employed to simulate more 
realistically the data management path from detector to memory, but the present 
profiles are sufficient to produce accurate data storage profiles. 
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“SOUTH ATLANTIC ANOMALY, 10 PROTONS/cm a /sec. >50 MeV 

NOTE: EACH 6 HOUR MISSION SEGMENT IS PRESENTED IN A THREE PAGE FORMAT. 


Figure 7. Option 1 timeline, TL-4. 
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F igur e 7 . (Co ntinued ) . 
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“SOUTH ATLANTIC ANOMALY, 10 PROTONS/cm 2 /sec, > 50 MeV 

NOTE: EACH 6 HOUR MISSION SEGMENT IS PRESENTED IN A THREE PAGE FORMAT. 


Figure 8. Option 2 timeline, TL-5. 
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E. Summary of Results and Statistical Analysis 

» 

A summary of quantitative results from the two DRM timelines is 
presented in Table 4. In addition to columns of statistics for options 1 and 2, 
a third column has been included to provide information relative to the prime 
viewing portion of the schedules to show how these statistics would change in 
the absence of secondary viewing. Recall that the prime viewing path is 
independent of the secondary path and is, therefore, the same in both options. 
Notice that although the slewing increases from 17 to 29 percent with the addi- 
tion of region C, the total lost time actually decreases from 46 to 43 percent 
because of the large drop in waiting time. This waiting time is idle time when 
the telescope is not observing due to one of the following reasons: (l) there is 
not enough time left at the end of the secondary interval to view another target, 
(2) viewing another target causes a drop in optimum path value, or (3) no more 
target opportunities exist in the interval. In efficient schedules, waiting time 
can be driven to very low values and will occur only at the end of secondary 
viewing intervals, before a slew from one secondary region to another or from 
a secondary region back to the primary. 

The foregoing statistics can be used to make a realistic estimate of the 
observation efficiency that can be achieved if baffling and other stray light pro- 
tection devices provide target intervals that span the entire sun side of each 
revolution. This improvement in target availability will allow the elimination of 
region C, so the slew and settle times would be similar to those of option 1. 

For the same reason, waiting time would be no more than that of option 2, and 
probably less. This would give an overall efficiency of ~73 percent, or, in 
other words, a 23 percent increase in efficiency due to the addition of secondary 
viewing. The efficiency of any given revolution, of course, depends on the 
amount of slewing required and duration of secondary target observation. 

Since no relation exists between the scientific worth and viewing time required 
on a target, however, efficiency is not as important as the number of targets 
completed and viewing should, therefore, not be biased with respect to time. 

Histograms of target magnitude, observation time, and slew time are 
presented in Figures 9 through 14. 
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TABLE 4. SUMMARY OF RESULTS FOR 2 DAYS (31 ORBITS) 

OF OBSERVATIONS 



Option 1 (TL-4) Option 2 (TL-5) Viewing Only 


Number of Targets Completed 


Average Number of Targets Completed 
per Revolution 


Average Number of Targets Viewed 
per Revolution 


Lost Time (%): 
Slew time 

Settle time 

Waiting time 


Observation Efficiency (%) 


Efficiency Range per Revolution (%) 


Number of Instrument Complements 
Scheduled 


Number of Slews: 

Intra-region A 
Intra-region 15 
Intra-region C 
A-B (B A) 
B-C 
C-A 


Average Slew Angle (deg): 

Intra-region A 
Intra-region B 
Intra-region C 
A-B (B-A) 

B-C 

C-A 

Overall 



48.0 58.7 


48.6 64.4 


46.5 50.4 
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TABLE 4. (Concluded) 


Statistic 

Option 1 (TL-4) 

Option 2(TL-5) 

Prime 

Viewing Only 

Average Slew Time (min): 




Intra-region A 

0.9 

0.9 

0.9 

Intra-region B 

1.5 

1.4 

- 

Intra-region C 

— 

1.1 

— 

A-B (B-A) 

7.6 

7.6 

— 

B-C 

-- 

10.8 

— 

C-A 

— 

8.7 

— 

Overall 

6.2 

7.2 

0.9 

Maximum Slew Angle (deg) 

94.1 

155.7 

2.4 

Minimum Slew Angie (deg) 

0.7 

0.7 

0.7 

Total Path Length (deg) 

5156.7 

9611.1 

10.0 

Average Slew Distance per Revolution (deg) 

166.3 

310.04 

0.32 

Unused Prime Viewing Time (%) 

- - ... 

0.51 

0.51 

0.5! 

Potential SAA Conflicts: 



i 

Percentage of Time 

1.63 

2.82 

i 

i 

' 1.35 

Number of Targets Affected 

10 

22 

2 

Number of Data Frames Taken 

103 

1 83 

mam 


IV. CONCLUSIONS 


Timeline studies in support of the LST, Spacelab and Integrated Mission 
Planning Activity (IMAP) programs have demonstrated that the SCREAM algo- 
rithm is a valuable analysis tool. It has been applied in parametric studies to 
measure the performance of viewing philosophies and slew rate and to generate 
high fidelity mission simulations. The results of these simulations have been 
successfully used by subsystem designers to construct pointing and control 
momentum histories and assess thermal control, electrical power, and data 
management systems. 
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VISUAL MAGNITUDE 


Figure 9. Distribution of targets scheduled in TL-4 according to magnitude. 

Timeline algorithms such as SCREAM have numerous applications and 
can be effectively used in all phases of a project as a focal point for system 
requirements and constraints, to assess the impacts of guideline and subsystem 
performance changes and to eliminate the confusion that often surrounds 
decisions whose implications are too complex to determine without realistic 
simulations. This in turn, could save much valuable time and effort in space- 
craft development. 
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VISUAL MAGNITUDE 


Figure 10. Distribution of targets scheduled in TL-5 
according to magnitude. 
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APPENDIX 


AN ALGORITHM FOR A SINGLE MACHINE SCHEDULING 
PROBLEM WITH SEQUENCE DEPENDENT SETUP 
TIMES AND SCHEDULING WINDOWS 


The selection of stellar targets for a telescope in Earth orbit gave rise 

to the following scheduling problem.* One machine is available for performing 

tasks between times T and T . The machine can perform only one task at 

s e 

a time. There exist n candidate tasks numbered 1 through n from which a 
schedule for the period is to be developed. The performance of task i requires 
W. time units and has a value of V. for i =* 1 > 2, . . . , n . In addition, task 

i l 

i is constrained to be performed within the time period that begins at time B. 


and ends at time . The machine must be adjusted before the beginning of 

each task. The time required to adjust the machine following completion of 
task i before task j begins as A„ . If task j is the first task in the 

sequence, then A is the time required to adjust the machine. The problem 

is to select a feasible sequence of tasks to perform between times T and T 
which will maximize 


Z 


l 


1=1 


x.v. 

1 1 


where X, = 1 if task i is selected and X. = 0 otherwise. 

1 i 

Let S , S , . . . , S represent a general sequence containing k 
d i d 

tasks where represents the number of the ith task in the sequence. Define 

to be the time when the performance of task d. is to begin. The problem 
i 

can now be stated formally as follows: 


*Guffin, O. T. and Haussler, J.B. ; Large Space Telescope Design Reference 
Mission Document. Marshall Space Flight Center, Sept. 1974 (Unpublished 
Report) . 
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From the n candidate tasks, select a sequence of tasks 
which maximizes 




n 


* - Z 

i=l 



subject to 



1. 

t = B for i = l, 
d. d. 

i i 

2* . • • > k . 

2. 

t , + W , = E , for 
d. d. d. 

ill 

i — 1, 2) ♦ * * 

3. 

t - T + A , . 

d ! s Od^ 


4. 

t + W + A d 
d. d. d. d+1 

i i i 

= t , for i 

d i+ i 

5. 

t , + W , ST . 
d, d, e 



This problem is very similar ,to the single machine job-shop problem 
with sequence dependent setup times [2] but it differs in two respects. First, 
the current problem is concerned with selecting an optimal subset of available 
tasks to schedule. In the job-shop problem, all available tasks are scheduled. 
The other distinction is that the performance of each task scheduled in the cur- 
rent problem is constrained to occur only within its particular time period. 

This constraint has no physical meaning in most job-shop problems. 

Since a job-shop problem can be formulated as a linear integer pro- 
gramming problem 1 3-6 1 , one would think that the current problem would have 
a similar formulation. However, consideration of the sequence dependent setup 
times and the selection of only a subset of the available tasks makes such a 
formulation very difficult. An enumeration algorithm similar to those used to 
solve zero-one integer problems was developed for solving the problem. 


Solution Procedure 

The algorithm implicitly enumerates only a portion of all feasible 
sequences while automatically discarding the remaining ones as nonpromising. 
The algorithm takes full advantage of the basic concepts employed by the Balas 
zero-one algorithm [7]. 
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The algorithm contains two distinct parts: 

1. A flexible enumeration scheme which records all sequences that have 
been considered and generates the remaining ones in a nonredundant fashion. 

2. A number of tests for excluding from further consideration those 
sequences which either do not have a better objective value than the best 
sequence found so far, or are not feasible, or both. 

Enumeration Scheme . An explanation of the enumeration scheme 
requires that the following definitions be introduced: 

( i) Partial Sequence (J) : This represents the first part of a sequence 
in which the first k tasks to be performed have been established. The notation 
J = (d t , d 2 , ... , d }, k < n is used to indicate that d. is the ith task to be 

K 1 

performed. A candidate task that is not contained in J is called a free task. 

(ii) Completion of J : This is a complete sequence whose first k terms 
are defined by the partial sequence J. The notation = {d t , d 2 , ... , d^, 

d. d } denotes a completion of J . J is complete in the sense that 

k+ 1 m e 

no free task can be added to the end of the sequence without violating a 

constraint. 

(iii) Fathomed Partial Sequence : A partial sequence is said to be 
fathomed if all of its completions can be discarded as nonpromising. Such 
completions may be either infeasible, or may not yield an improved value of 
the objective function, or both. Tests are applied to determine when a partial 
sequence is fathomed. These tests are called fathoming tests. 

The enumeration scheme is based on the characteristic which distin- 
guishes between two distinct sequences. Suppose G={d t , d 2 , ... , d } and 

K. 

H = {h t , h z , ... , h^} are two sequences. Then G = H if and only if m = k 
and d. - h. for i - 1 , 2 , ... , k. The enumeration scheme considers all 

ii 

possible sequences by generating a nonredundant sequence of partial sequences 
J 0 , J t , J 2 , ... , J^, **• where differs from each J 0 , Jj, ... , 

in at least one position. 
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A number of variables will now be defined for use in explaining the 
procedure for generating a nonredundant sequence of partial sequences. 

the number' of positions in the current partial sequence J 

K 

= {d, , dg , . . . , d^} is the current partial sequence 

number of free tasks which have been found to cause infeasibility 
when assigned to position i of J 

K 

jth free task which cannot be assigned to a position for feasibility 
reasons; Q defines the specific tasks which are associated with 
P. 

l 

number of candidate tasks that have been considered for ith posi- 
tion in Jj since the last time a task change occurred in any of 

positions 1, 2, . . . , i - 1 

ith candidate task which has been considered for jth position in 
Jj_ since a task change occurred in any of positions 1, 2, ... , 

j - 1 ; variables It and H record which sequences have already 
been considered 

= { 1, 2, . . . , n} is the set of all candidate tasks 
The steps in the enumeration scheme. can now be listed. 

Step 1; Initiation 

J k ~ <p , r = 0 , k = 0 and P = it. = 0 for ieN . 

Proceed to step 2. 

Step 2: Apply fathoming test. Can J, be fathomed? If no, continue 

k 

to step 3. If yes and r > 0, go to step 5. Otherwise, go to 
step 8. 

Step 3: Forward Move : Increment r by one. If r = 1 , set = 0 ; 

otherwise set P = P . Define C to be the set of all 
r r- 1 

free tasks excluding those defined by for i = 1, 2, . . . , 
P r « If C is empty, go to step 6. Otherwise, for each 


Q. 


R. 


H 


ij 


N 
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a eC , determine whether the sequence consisting of J 

k 

with a augmented as its last element is feasible* For each 
a which causes the augmented sequence to be infeasible, 
increment P y by one and set Q p =a. Is there an a€ C 

r 

which retains feasibility in augmented sequence? If yes, go 
to step 4. If no, go to step 6. 

Step 4: R r = 1 and = a . Form J^ +1 by augmenting a to J 

k 

as its last element on the right. Add one to k and go back to 
step 2. 

St e P 5: Lateral Move : Define C to be the set of all free tasks 

excluding those defined by Q. for i= 1, 2 P and 

1 r 

te r ter j - 1 , 2 , ... , R^ , If C is empty, go to step 6. 

Otherwise, for each o<€C, determine whether the feasibility 
of would be maintained if d^ is replaced by a . For 

each a that would cause infeasibility, increment P by one 

r 

and set Q p = a . Is there an neC which would maintain 
r 

feasibility? If yes, go to step 7. If no, go to step 6. 

Step 6: Backward Move : Decrement r by one. This deletes the 

right-most position of J^. Is r > 0 ? If yes, go back to 

step 5 and attempt to select a replacement for d . If no, 
proceed to step 8. r 

Step 7: Form J fc+1 by replacing d f with a from step 5. Increment 

R r by one and set j = R f . Set H « a . This records the 

consideration of task a for the rth position in current 
partial sequence. Increment k by one and return to step 2. 

Step 8: Stop. All possible sequences have now been implicitly 

enumerated. 
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Note that the scheme is characterized by three basic steps: The 
forwai-d move augments a free task on the right of the current partial solution; 
the lateral move replaces the task in the right-most. position; and the backward 
move drops the right-most position. Schemes of this type are referred to as 
backtracking procedures. 

The definitions of the variables P and Q may leave some doubt as to 
how they record infeasibility. Q. for j = 1 , 2, P defines the tasks which 

cannot feasibly occupy position i . Q. for j = P. + 1 , P. + 2 , ... , P. 

j 11 H* 1 

defines the tasks which could feasibly occupy position i but would cause 

infeasibility in position i + 1 . Notice that any task which would be infeasible 

in the ith position would also be infeasible in positions i + 1 , i + 2 , .... That 

is, if a task cannot be performed between times t, + W and T , then it 

d. d. e 

l l 

certainly cannot be performed between times T , + W and T where 

d i+ m d i + m 

m is a positive integer. This is the reason that P = P , at initiation of 

r r— 1 

step 3. 


The validity of the enumeration scheme is ascertained by the following 
theorem: 

Theorem : (i) The completions of all future partial sequences do not 
duplicate any of the completions which have already been fathomed. 

(ii) The scheme terminates only after all feasible sequences 

are considered. 

Before providing a proof, a general description of the way in which the 
scheme operates will be given. In general, the scheme considers a group of 
feasible sequences at one time as opposed to considering one sequence at a 
time. A group of sequences is defined by a partial sequence in which the first 
r positions have already been assigned specific tasks. The completions of a 
partial sequence define the group of sequences to be considered. If the com- 
pletions of a partial sequence are not fathomed, the number of completions 
under consideration is reduced by augmenting one more task to the partial 
sequence. The augmentation procedure is continued until the number of com- 
pletions is small enough to be fathomed. At this point, the task in the last 
(rth) position of the partial sequence is replaced if a feasible and nonredundant 
replacement can be found. A replacement is nonredundant if it has not pre- 
viously occupied this position since the task in the (r - l) position has been 
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replaced. The mechanics, of ithe ■ scheme in^u^e ..fhat J& e , tasks , in. positions 
1 , 2, ... , r - 2 have not been replaced if the task in position r - 1 has 


if a feasible and nonredundant replacement is available. If a replacement is 
available, the new set of completions being considered is nonredundant. Other- 
wise, the rth position is deopped and r decremented by one. The dropping of 
the last position is continued until either a position is reached for which a 
feasible and nonredundant replacement is available or r reaches the value of 
zero. When r becomes zero, all feasible completions have been implicitly 
enumerated. 

Proof of Theorem : The proof will be by induction. The eight steps of 
the enumeration scheme will be referred to during the proof. The theorem is 
obviously true for n = 1. Assuming that it is true for n = u , we proceed to 
show chat it is true for n = u + 1 . 

If the initial null partial sequence can be fathomed at step 2 , the enumera- 
tion shceme is not used. So, we assume the process proceeds to step 3 after 
step 2. At step 3, a test is performed on each candidate task to determine if it 
can be feasibly performed as the first task in a sequence. Variables P and Q 
record which task cannot feasibly occupy the first position. If > 0, the 
problem reduces to one containing, at most, u tasks. This follows from the 
definition of the set C in steps 2 and 5. C can contain no more than u + 1 - P 

r 

tasks where ^ Pj for r > 0. Thus, if the method can fail, it can only 

do so when P t = 0 at the end of step 3. In this case any one of the u + 1 tasks 
could feasibly occupy the first position. One is selected and assigned to the 
first position to form Jj at step 4. The assignment of this particular task to 
the first position is recorded by variables R and H . The problem is now 
temporarily reduced to one containing u tasks since the task just assigned is 
removed from the list of free tasks. 

By the induction hypothesis, the enumeration scheme considers all 
feasible sequences of the u task problem without duplication. Termination of 
this problem occurs at step 6 when r is reduced once again to a value of 1. 

The process now proceeds to step 5 to locate a replacement for the task initially 
assigned to the first position. 


not been replaced. Thus replacing the task in the "rth position' insures that 1 
the next completions Ipeing p6nsid,ered are'not redundant; ( If a feasible arid 
nonredundant replacement for the ' rtti* Vositiob cannot be found," then all feasi 
ble completions of the, current . partial .sequence with l its last position deleted 
have been fathomed. The,’ last position of curreht pa’r tial se q ue'ric e is' dropped 
and r is decremented by one. The task in the new rth position is replaced 
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At step 5, the candidate tasks from which the replacement is to be 
chosen are defined by the set C. Since we established at step 3 that Pj = 0, 

C contains all u + 1 tasks except the one recorded by the R and H variables. 
(Note that will retain a value of zero because we assumed above that each 
of the u + 1 tasks could feasibly occupy the first position in a sequence. ) This 
means that the task chosen for the first sequence position will be different for 
the one initially assigned to this position. Therefore, all completions of J t 
which have already been fathomed will not be duplicated by any completions of 
the new partial sequence that is created at step 7. 

The assignment of another task to the first position has once again 
temporarily reduced the problem to one involving u tasks. The process now 
returns to step 2 and begins consideration of this u task problem. By hypoth- 
esis, the process considers all feasible sequences of this u task problem 
without duplication. Termination occurs at step G when r is reduced to a 
value of 1. Once again a replacement for the task in the first position is 
selected at step 5. Since the first two variables considered for the first posi- 
tion were recorded by the R and H variables at steps 4 and 7, the task 
chosen at step 5 will be one not previously considered for the first position. 
Therefore, the completions of the new partial sequence will not duplicate any 
of those that have already been fathomed. The problem has again reduced to 
one containing u tasks at step 7. 

The scheme has now begun to follow a set iterative pattern. The ter- 
mination of a u task problem occurs at step 6 when r is decremented to a 
value of 1. A replacement for the task assigned to the first position is chosen 
at step 5. At step 7, the chosen task is assigned the first position and its 
assignment is recorded by R and H variables. This establishes a new u 
task problem. 

To establish that the scheme considers all feasible sequences in a non- 
redundant way before terminating, notice that at step 5 the replacement task is 
always chosen from among those tasks not previously assigned to the first 
position. This is assured by the use of the R and H variables in establishing 
the set C. Since P t = 0 , C always contains all tasks not previously assigned 
to the first position. This means that C will be empty and cause the scheme to 
terminate only after all u + 1 tasks have been previously assigned to the first 
position. This completes the proof, since the first sequence position can only 
be assigned in the u + 1 distinct ways that have been considered in this process, 
and, by assumption, in each of the u + 1 cases that the method handles the 
resulting u task problem as specified by the theorem. 
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Fathoming Tests. The fathoming tests determine whether J has 

any promising completions. If none exist, the lateral move of the enumeration 
scheme occurs, and additional fathoming tests determine which task should be 
considered next for the rth position in . If all completions of cannot 

be shown to be nonpromising, other fathoming tests determine which task should 
be augmented to during the forward move . 


Before proceeding to describe the fathoming test applied at step 2, 

additional variables need to be defined. Let Z, be the value of the best 

b , 

r 

sequence at present. Define Z^ to be the value of J , that is, Z^ = ^ 

i=l 


three 



Define T to be time when performance of task d^ is completed. These 
three variables will now be used to explain the fathoming mechanism. 


As stated above, the fathoming test determines whether J has any 

K 

promising completions. The determination is made by computing the highest 

possible value, Z, , any oGiiiiJieiiuii of J. couid attain, it Z f ^ Z_ , then J. 

^ k k kb k 

will be fathomed. Otherwise, J, is not fathomed. 

k 

is computed by first determining which free tasks would be feasible 
candidates for occupying positions r + 1, r + 2, ... in any completion of J . 

K 

Since any task that couid feasibly occupy position r + 1 could possibly occupy 
positions r + 2 , r + 3 , ... , the list of candidates is restricted to those which 
couid occupy position r + l . Next, these candidates are ordered according to 
their ratio of completion value, , to performance time, W. . The candidate 

having the largest ratio is placed in first position of order. For simplicity of 

notation, let 1, 2, 3 m be the ordered list of candidate tasks. Next, 

compute M = min {A ^ . |i = 1, 2, ... , m) . M is the minimum setup time 

r 

between task d^ and any one candidate task. Consider the sequence which 

places task 1 in position r + 1 of J with a setup time of M, task 2 in 

position r + 2 , etc. , until all of the time between T and T has been filled. 

re 

Define U = T - T . Then tasks 1, 2, , are placed respectively in 

positions r + 1 , r + 2 , ... until X - U - M - W t - W 2 - . . . - W ^ 0 and 
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X - W <0. Task q + 1 is placed in position r + q + 1 and a value of 
q+1 

X( V /W ) is assumed for its performance value. The value of the best 
v q+1 q+T ' 

possible completion of J could not exceed the value of this sequence. 

Z, is computed as follows: 
k 


Z, = Z, + V V. + (X/W , J V 
k k i q+1 q+1 


For example, suppose M = 1, U = 17, W t - 7, W 2 - 7, W 3 =5, W 4 = 6, V t = 24, 
y 2 = 14 , v 3 = 8, and V 4 = 6 . Then tasks 1, 2 and 3 are placed in positions 
r+ 1 , r + 2, and r + 3 , respectively. Since there is only time to perform 
task 3 for two time units, its performance value is reduced to (2/5) 8 = 3.2. 
Then 


Z. = Z + 24 + 14 +3.2 
k k 


The execution of either the forward move (step 3) or the lateral move 
( step 5) requires that a task be chosen for augmentation or replacement, 
respectively. In either case, a heuristic rule for selecting the most 
’’promising” task is used. Intuitively, this task should provide a maximum 
amount of value and require a minimum amount of setup, waiting, and per- 
formance time. This is equivalent to choosing the task with the maximum ratio 
of value (V.) to total (setup plus waiting plus performance) time. If the period 

of time between T „ and the time the ith candidate task can begin (B.) is 
r-1 i 

longer then the setup time (A^ ^ .) , the machine would be idle for B. - A f ^ . 

T time units. The waiting time associated with the ith candidate task is 
r-1 

therefore max( 0, B. - A - T^ ^) . The ratio for the ith candidate task 
is 


V 

l 

w + A . + max( 0, B. - A ~ . - T H 
i r-1 l i r-1 l r-l 
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An Example. The operation of a machine is to be scheduled between 

times T = 0 and~ T =20. There exist three candidate tasks from which a 
s e 

schedule is to be developed. The following matrix lists the characteristics of 
each task: 


Task 

W. 

i 

V 

i 

B. 

l 

E. 

l 

A 0j 

1 

8 

12 

11 

22 

3 

2 

5 

3 

3 

11 

4 

3 

7 

5 

4 

17 

2 


The following matrix defines the sequence dependent setup time (A..): 

The eight steps of the enumeration schemes 
will be referred to as the algorithm is applied. 


Task j 


c i 

1 

4. 

<rt 

0 

i 

0 

4 

2 

2 

3 

0 

1 

3 



1 

2 

j 

0 


Step 1: 


Initiate variables. 

Set r = 0, k = 0, J = <P , 
Z 0 = 0 , z b = , z k = 0 

and T = T =0. Proceed 
r s 

to step 2. 


Step 2: Apply fathoming test. Since Z 0 > Z^, J 0 cannot be fathomed. 

Continue to step a. 

Step 3: Select a task to augment to J 0 . Set r = 1 , Pj = 0 and C = 

{ 1, 2, 3}. For each yeC , determine whether the sequence 

consisting of J 0 with y as its only task is feasible. The 

sequence consisting of task 1 is feasible since A 01 + T 0 < 

and < T < Ej . Also, the sequence consisting of 

c 

task 2 is feasible since T 0 + > B 2 and T 0 + A 02 + W 2 < 

E 2 < T , Similarly, the sequence consisting of task 3 is 
O 

feasible. 

The value to total time ratio 
V. 

l 

W. + A , . + max( 0, B. - A , . - T T 
i r- 1 i i r-li r- r 
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Step 4: 
Step 2: 

Step 5: 

Step 7: 
Step 2: 

Step 3: 


is now computed for each yeC to determine which of the 
three tasks to augment to J 0 . The ratios for tasks 1, 2, and 
3 are 


12 

8+3+8 


3 

5+4+0 


and 


5 

7 + 2+2 


respectively. 


Task 1 is selected since it has the largest ratio. Proceed to 
step 4. 


Augment task 1 to J 0 . Set k - 1, R t = 1, = 1, = 

{ 1} , Z t = 12, t t = 11, and T t = 19 . Go to step 2. 


Apply fathoming test. First determine the free tasks that 
could be feasibly augmented to Jj . Tasks 2 and 3 are the 
free tasks, but neither of them could be feasibly augmented 
to Jj because T t > E 2 and Ti > E 3 . Therefore, J t is a 
complete sequence. Set Z^=Z 1 =12. J t is fathomed since 

it has no more promising completions. Proceed to step 5. 


Select a replacement for task 1 in J ( . Since Pj = 0, R t = 1 
and H lt = 1, C = { 2, 3} . Tasks 2 and 3 are both feasible 
replacements. The value to total time ratios for tasks 2 and 
3 are 3/9 and 5/ll, respectively. Therefore, task 3 is 
chosen to replace task 1. Proceed to step 7. 

Set k = 2, J 2 = » Rj = 2, n 2 ^ — 3, t^ — 4, Tj - 11, and 

Z 2 = 5. Go to step 2. 


Apply fathoming test. First determine the free tasks that 
could be feasibly augmented to J 2 . Tasks 1 and 2 are the 
free tasks. Task 1 can be feasibly augmented since T t + 

A 31 > Bj and Tj + A 3i + < F t . Task 2 cannot be 

feasibly augmented since Tj + A 32 > _E 2 . Compute Z 2 = Z 2 + 

V t = 17 . Jo is not fathomed since Z, > Z, . Continue to 
. „ 1 b 

step 3. 

Set r = 2, P 2 = Pj - 0, and C = { 1, 2} . In the previous step, 
we learned that task 1 could be feasibly augmented to J 2 , 
but task 2 could not be feasibly augmented. Set P 2 = 1 and 
= 2 . Proceed to step 4. 
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Step 4: 


v 


Set k = 3, J 3 = {3, 1} , Rg = 1, H 12 = 1, Z 3 = 17, tj = 12, and 
T 2 = 20. Go to step 2. 

Step 2: Task 2 is the only free task and it cannot be feasibly augmented 

to J 3 since T 2 > E 2 . This means that J 3 is a complete 

sequence, so it is fathomed. Set = Z 3 = 17 and go to 
step 5. 

Step 5: Since P 2 = 1, = 2, Rg = 1, and H 12 = 1» C is null. 

Continue to step 6. 

Step 6: Set r= 1, J 3 = {3}, t 3 = 4, T t = 11, and to to step 5. 

Step 5: Attempt to find a replacement for task 3 in J 3 . Since Rj = 

2, H lt = 1, H 21 = 3, and P t = 0, C = {2} . Task 2 can 
feasibly replace task 3. Since it is the only task in C , it is 
selected to replace task 3. Proceed to step 7. 

Setp 7: Set k = 4, J 4 = {2>, R< = 3. H S1 = 2, t 2 = 4, T t = 0, and 

Z 4 = 3. Go to step 2. 

Step 2 : Apply fathoming test. First determine the free tasks that 

could be feasibly augmented to J 4 . Tasks 1 and 3 both 
satisfy this condition. Next, order these two tasks according 
their V./W. ratios. Task 1 is placed in first position since 

it has the largest ratio. The minimum setup time between 
task 2 and either task 1 or task 3 is 1, so M - 1. Compute 
U = T - Tj = 11. Since X = U - M - W t > 0 and X - W 3 < 0 , 

Z 4 = Z 4 + V 1 + X(V 3 )/W 3 = 16 + 3/7 


J 4 is fathomed because Z 4 < = 17. Go to step 5. 

Step 5: Attempt to find a replacement for task 2 in J. . C is null 

because Rj = 3, H 14 = 1, H 21 = 3, and H 3l = 2. Continue to 
step 6. 

Step 6: Set r = 0. Proceed to step 8. 

Step 8: Stop. All possible sequences have now been implicitly 

enumerated. 
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The optimum sequence is given by J 3 = { 3, l} with = 17. Notice 
that it was identified without enumerating ail feasible sequences. 


Computational Experience 

The algorithm is currently being used to support the design of a tele- 
scope to be operated in Earth orbit. Typical observational timelines are 
developed and used in the design of the telescope' s pointing, maneuvering, 
data management, and other subsystems. 

An observational timeline specifies the time sequence in which targets 
are to be viewed over several orbits of the telescope. The telescope is shaded 
from sunlight by the Earth for a portion of each orbit. This period is always 
used to observe a very faint target. Each faint target scheduled must be viewed 
during the shadow portion of several consecutive orbits to complete data col- 
lection. Brighter targets are viewed during the remainder Of each orbit. 

Since the brighter targets require only a short time span to complete data 
collection, several of them may be viewed during the sunlight portion of an 
orbit. The problem of selecting a sequence of bright tai'gets to view during the 
sunlight portion of one orbit coincides with the sequencing problem presented in 
this paper. The scheduling routine used to generate a timeline for several 
orbits employs the algorithm presented herein to develop a timeline for the 
sunlight portion of each orbit. 

The algorithm was implemented in a Fortran IV code for a Univac 1108 
computer. To obtain computational experience, the code was used to generate 
a timeline for the sunlight portion of 10 typical orbits of the telescope. Each of 
the 10 timelines was generated from a different list of 200 bright targets. The 
number of targets per optimum timeline ranged between 6 and 10. The com- 
putational time, excluding input and output times, per timeline ranged from 
0.61 sec to 0.80 sec with a mean of 0. 71 sec. The same 10 problems were 
rerun after the size of each candidate list of bright targets had been reduced 
to 150. The computational times ranged from 0. 38 to 0. 42 sec with a mean of 
0.40 sec. Just as one would expect, the size of the target list has a significant 
influence on computational time. Additional computational experience indicates 
that computational time also increases as thenumber of targets in the optimal 
sequence increases. 
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Other Potential Applications 

The algorithm is presently being considered for use in scheduling the 
daily activities of a crewman during a space mission. There usually exist 
many more daily tasks, not all of which are mandatory, than a crewman can 
perform. The value or priority of the different tasks vary. In addition, the 
performances of many tasks are constrained to occur within particular time 
windows due to target availability, housekeeping requirements, etc. (If task i 
does not have a scheduling window, set B. = T and E. = T . ) Frequently, as 

L S It? 

during the Skylab mission, a crewman will have to refer to a manual (set him- 
self up) before performing the task. These constraints are considered by the 
algorithm. However, this problem usually has additional constraints such as 
equipment conflicts, consumable constraints, and crewman availability. An 
extension of the algorithm is presently being developed to consider these con- 
straints. Similarly, one person scheduling problems, such as the sharing of 
a quality inspector among several departments, may exist outside of NASA. 
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